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1. Introduction 

Direct searches and electroweak precision tests suggest that a Standard Model Higgs boson 
is light, and that it should decay predominantly into a bottom quark (bb) pair. Inclusive 
searches at the LHC of a Higgs resonance in the bottom-pair invariant mass distribution 
are unfortunately hampered by a very large QCD background. Direct searches are focusing 
on other decay channels, such as H — > 77, with a much smaller branching ratio. In view 
of this, perturbative corrections to the H — >■ bb decay are most interesting for the inclusive 
rate, due to its contribution to the total decay width and the branching ratios of other 
decays. A remarkable effort has been made for a precise calculation of T H ^ b i and now 
QCD corrections are known to 0(a A s ) [1-9]. 

Recently, confidence has grown that the H — > bb decay can be observed at the LHC 
directly in events where the Higgs boson is produced in association with a massive vector 
boson (W, Z) [10] or a tt pair [11]. Backgrounds from ti,V+jj and multi-jet production are 
challenging. However, the excellent 6— jet tagging of the ATLAS and CMS detectors as well 
as sophisticated selections of jet events described in [10,11], render these channels hopeful, 
as explicitly demonstrated by the ATLAS collaboration in a full detector simulation analysis 
for the ZH case [12]. For this channel, further suppression of the background can be 
obtained with the sub-jet algorithms of [13]. 

These search strategies rely on a selection of phase space corners which are rich in 
potential Higgs events. An accurate modeling of QCD radiation is necessary in order 
to assess the efficiency of these selections. This motivates the computation of the fully 
differential H — > bb decay rate at next-to-next-to-leading-order (NNLO) in perturbative 
QCD. 

Next-to-leading-order (NLO) computations [14-26] are currently performed with very 
well automated methods [27-35]. Obtaining fully differential cross-sections and decay rates 
at one order higher in the perturbative expansion requires the solution of new challenging 
problems. In the last decade rapid progress [36-60] with important results [61-73] has 
been achieved. However, further developments of methods and new ideas are necessary for 
efficient cancelations of infrared singularities and evaluations of novel two-loop amplitudes 
in more complicated LHC processes. 

In Ref. [74] we introduced a new method to extract the singularities of multi-dimensional 
loop and phase space integrals which emerge at NNLO. The method employs non-linear 
transformations for the factorization of overlapping singularities which in turn allow for 
an efficient numerical evaluation. In this article, we demonstrate the method in a realistic 
physical example. We compute for the first time the fully differential decay rate for H — > bb 
at NNLO. 

The paper is organized as follows. In Section 2 we fix our notation and set up the 
calculation. In Section 3 we present the matrix elements for the processes which contribute 
to the decay width through NNLO. In Section 4 we describe the parameterization of the 
phase space for the final state partons. In Section 5 and Section 6 we explain how our 
method leads to a simple factorizable form for the singularities in the integrals which 
emerge at NNLO. In Section 7 we present the inclusive width and examples of differential 
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observables with our numerical code. 



2. Notation 

We compute the partial width ^H^-blx [J\ through NNLO in perturbative QCD for any- 
infrared safe observable J, such as an appropriate jet- algorithm, in the decay of a Higgs 
boson to a pair of bottom-quarks, 



H -> b + b + X. 



(2.1) 



We consider that the Higgs boson couples to a bottom-quark pair with an unrenormalized 
Yukawa coupling y? . Due to the largely different mass of the Higgs boson, m#, and the 
small bottom-quark, m&, we treat the latter as independent of the Yukawa coupling and set 
it to zero, m& = 0, in matrix elements and phase space integrals. In addition, we neglect 
top-quarks both in virtual corrections and in the renormalization procedure and consider 
rif = 5 light quark flavors. For the renormalization of ultraviolet (UV) divergences we 
work in the MS-scheme where the bare couplings, yf, and gf = Airaf , are related to 
their renormalized counterparts through the following relations: 
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(2.2) 
(2.3) 



The renormalized strong coupling and Higgs Yukawa coupling depend implicitly on the 
renormalization scale /x, 

a s = a s (n) yb = yb{n). (2.4) 



The SU(N) Casimirs are given by 



C A = N, C f = 



N 2 



2N 



T F = 



(2.5) 



where TV = 3 is the number of quark colors. 

Through NNLO, the partial decay width receives contributions from the following 
partonic processes: 

• H(ph) — >• b(pi)b(p2) at two loops, 

• H(ph) — > b{p\)b{p2)g{pz) at one loop and 

• H(p H ) -> b(p 1 ^p2)g(p3)9(p4),b{pi)b(p2)b(p3)b(p4),b(p 1 )b(p 2 )q(p3)q(p4) with q ^ b 
at tree level, 
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where in parenthesis we indicate the four-momenta of the corresponding particles. The 
width is the sum 
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In the above, we integrate the appropriate matrix elements computed in conventional 
dimensional regularization (CDR) over the phase space in D = 4 — 2e dimensions of the 
final-state partons and the function J defining an infrared safe observable. The phase 
space measure reads 



n \ n , D 



d* n ^ {2k)°SW [ph-^pAU 

V i= i / i=i ( 27r ) 



Lorentz invariants are defined as 



mjj — pjj, Sij — 2pi.pj, Sijk — Sij + Sjfc + Sjk. 



(2.8) 



(2.9) 



Regarding the analytic continuations of complex functions presented in this paper the 
positive "epsilon prescription" must be used, since all occurring mass scales are external, i.e. 



m 2 H — > m 2 H + ie, Sjj — >• + ie. 



(2.10) 



where e is a small positive parameter. 



3. Matrix elements 

In this section we present the square of the scattering amplitudes needed for the compu- 
tation of the H — > bb decay width through 0(a s (fj,) 2 ). These can be easily derived with 
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modern methods for computing Feynman diagrams and have been ingredients of many 
prior calculations. For example, the same matrix elements enter the calculation of the 
NNLO total cross-section bb — > H [75]. We have made an independent computation and 
present them here for completeness. For the generation of the matrix elements we used 
QGRAF [76]. For further symbolic manipulations, such as color and Dirac algebra, we 
used FORM [77] and MAPLE [78]. We reduced one and two-loop amplitudes to master 
integrals using the method of integration by parts [79, 80] and the Laporta algorithm [81] 
with AIR [82]. 

3.1 Decay to two partons 

For the process 

H( PH ) -> 6(pi) + b(p 2 ) 
we require up to 2-loop corrections: 
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and the two-loop amplitude interfered with the tree amplitude 
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with ^ = ln(4). 

3.2 Decay to three part 

For the process 

H(p H ) -> b(pi) + b(p 2 ) + g(p 3 ) 
we require up to 1-loop corrections: 

«sQ) .(0) / 2 ,(1) 
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is the quark gluon splitting kernel. The one-loop correction can be expressed as 
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There are only two different master integrals which appear here. The bubble 
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3.3 Decay to four partons 

For the decay to four partons we need 

H(ph) -»■ 6(pi) + b{p 2 ) + i(p 3 ) + j(Pi), (u) € {(gg), (66), (55)} 

at tree-level 
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The H — > 66gg amplitude can be expressed as 
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The form factor B(pi,p2,P3,P4) is presented in Appendix A. The squared amplitude 
corresponding to the Higgs decaying into a bb plus two extra gluons can be expressed as 

■*£Uw = 16?r4 (^r) ^ {° fA ™ 99 + N C 2 F B Hbhg + C F (l + 2NC F )C mlgg ) (3.20) 
where 

A - = A {0) + e4 {1) + e 2 4 {2) 

Hbbgg Hbbgg ~ Hbbgg ~ Hbbgg 

B - - R (0) + f f?W + f 2o(2) 
Hbbgg Hbbgg ~ l ~ Hbbgg ~ T ~ Hbbgg 

C Hbbgg = C Hbbgg + eC Hbbgg- ( 3 ' 21 ) 

The form factors A ( ^ l7 , MlL and C^'L can be found in Appendix A. 

Hbbgg Hbbgg Hbbgg 1 1 

4. Integration over phase space 

The matrix elements of Section 3 are integrated over the phase space of the final-state 
partons. For convenience, we perform our calculation in the rest frame of the Higgs boson. 
Our results can be trivially extended to any frame of reference. The matrix-element for 
the 1—7-2 process H — > bb is independent of any phase space integration variables, and it 
is multiplied by the phase space volume: 

$2 "j d$2 "(2^^ (%) "Svr r(2-2e) " (4 ' 1} 

The matrix elements for the 1 — > 3 processes are integrated over the corresponding 
phase space with volume: 



*3 = (i^7# - e) X dX^X^rX^. (4.2) 



The invariants take the simple form: 

si2 = m H X 2 

S13 = VUhMXx 

«23 = rn 2 H X 2 Xi (4.3) 
In the above, we have introduced the shorthand notation 

A = 1 — A. (4.4) 
We parameterize the phase space for the 1 — >■ 4 processes as follows: 

$4 = ^4 ! dX^dX^dXdX^X^^^X^XiX^ [sin(A 5 vr)]~ 2e (4.5) 
J o 

where 

iV 4 = (^) 2 ' 3£2 ' 13+8£ f 4+3e(1 ~ 2e) . (4.6) 
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The invariants in this parameterization are 

S234 = rn 2 H Xi 
S34 = m 2 H X 1 X 2 
«23 = m^AiA 2 A 4 
S24 = 171^X1X2X4 
S12 = m^AiA 2 A 3 



and 



S134 = m#(A 2 + A3A1A2) 



■513 = mjjXi A4A3 + A2A3A4 + 2 cos(A57r)\/ A2A3A3A4A4 



s 14 = m 2 H Xi A3A4 + A2A3A4 - 2 cos(A57r)\/ A2A3A3A4A4 



(4.7) 



(4.8) 



This parameterization is very similar to the ones used in Refs [37,57]. 

In the following sections we shall describe the method of non-linear mappings that we 
have employed in order to perform the non-trivial phase space integrations over the matrix 
elements at NNLO. Specifically, we shall discuss the "real-virtual" contributions which 
require both an integration over the three-parton phase space and an one-loop amplitude, 
and the "double-real" contributions which require the integration of tree matrix elements 
over the four-parton phase space. 



5. Real- Virtual NNLO contribution 

In the one loop amplitude for H — > bbg we will face integrals like 
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to extract the singularities. To tackle this difficulty we use Euler's integral representation of 
the hypergeometric function, and apply a non-linear mapping to disentangle an overlapping 
singularity of the integration variable in the hypergeometric representation and the phase 
space variables which generate 513,523. Defining 

^ 2^1(1, -e, 1 - e, -z) 
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we apply the transformation, 

x 3 ^ (3(x 3 ,l,z), P(x,A,B):= 
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which yields 

F(z) = (1 + zf £ dx 3 x^ (l - € . (5.5) 

This representation of the hypergeometric function has finite values at the points where 
we need to make a subtraction, 

2 F 1 (-e,-e,l-e,0) = 1 

2 F 1 (-e, -e, 1 - e, 1) = T(l + e)r(l - e), (5.6) 

and we can obtain the Laurent expansion in e of the real-virtual contribution easily with 
a "plus-prescription" subtraction. For non-special values of the kinematic variables the 
e— expansion of hypergeometric function reads: 



iFi (-e, -e, 1 - e,z) = 



The mapping of Eq. 5.4 simply re-derives a well known identity 

2 F 1 (a,b,c;z) = (1 - z^^F^c - a, b, c; ^—^ ), (5.8) 

and there would be no need of it had we started with the following representation of the 
box master integral: 

„ , 2r(l + e)r(l-e) 2 / . , , 

Box(u,t,M 2 ) = e2r(1 _ 2e) ("*) 

- {-M 2 )-%ut + sM^SF.i-e, -e, 1 - e, * ) (5.9) 

ut + sM z 

+ (-« - s) e 2 i ? i(-e, -e, 1 - e, — ^— ) + (-/ - s) e 2 i ? i(-e, -e, 1 - e, -^-) 

w + s i + s 

On the other hand, more complicated master integrals for one-loop amplitudes do not 
always have representations in terms of hypergeometric functions with known expansions in 
e and variable limits. However, it is always possible to derive a Feynman representation or a 
hypergeometric representation (via a Mellin-Barnes representation) for them and attempt 
a direct numerical evaluation of the combined phase space and virtual multiple integral 
using non-linear transformations to factorize the singularities. We have compared the two 
approaches: 

1. Evaluate numerically the coefficients of the Laurent expansion for the two-dimensional 
phase space integral after we express the loop amplitude and its hypergeometric func- 
tions in terms of polylogarithms, 

2. Evaluate numerically the coefficients of the Laurent expansion for the three-dimensional 
combined phase space and loop integral. 

We find that the three-dimensional integration is as fast as the two-dimensional integration 
where we require the evaluation of polylogarithmic functions. 
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6. Integration of double-real contributions 



phase space integrals for four-parton processes develop singularities at multiple soft and 
collinear kinematic configurations. These singularities are overlapping. In this section, we 
demonstrate how they can be factorized with non-linear transformations. We begin the 
discussion with an integral with four singular denominators, 

h [J] = [ d* 4 J{Pl ' P2 ^ P4) . (6.1) 

where the function J represents a non singular function of the partonic momenta, composed 
of numerators in the squared matrix elements and the infrared-safe observable J . We 
disentangle the overlapping singularities using 

A 3 ^a(A 3 ,A 4 ,l), (6.2) 

where 

a( ^ B): =iiW < 6 - 3 » 

We note that x i-> a(x, A, B) 4>iH a(x, B, A). After this transformation, the expansion 
in e is straightforward with simple subtractions. For the special case J = 1, we obtain 



h [J = l] = N 4 



m 7 «m + 8 ' 554 ( 2 ) °- 171 °( 2 ) . 0-34657(2) , °- 25 
61.76(2) + - 2 + -3 + — r 



(6.4) 



Other integrals may be mapped to I\ of Eq. 6.1 with a relabeling of the partonic 
momenta. For example, the integral 

h = I dg /fo'**'**^) , (6.5) 

J Si4S23'Sl34'S234 

which has an apparently more complicated singularity structure due to S14 in the denomi- 
nator, is transformed to the integral of Eq. 6.1 by exchanging P2 and P4. 
A second class of integrals has the following denominator structure: 

I 2 [J] = [ d * A J (Pi>P*P*>P*) ( 6 . 6 ) 
J S13S23S134S234 

Notice, that the integrand is singular not only at the boundaries of the integration variables 
\i, but it also develops a "line singularity" inside the integration volume. We eliminate 
this by casting the integral in the form 

T r Tl [ « (J(Pl,P2,P3,P4:) + J(P2,P1,P3,P4))S24: (R 7 , 

h[J\ = / ; > • (6.7) 

J S23Sl34S234(Sl3S24 + 514^23) 

It is easy to see that Eq. 6.7 is equivalent to hlJ], by exchanging the momenta p 1 o p 2 i n 
the term with J(p2,Pi,P3,P4:)- The denominator 



S13S24+S14S23 = m|^AiAiA2 2A3A4A4 + A2A3(A| + A 4 ) + 2cos(7rA5)(l - 2X4) \/ A2A3A3A4A4 

(6.8) 
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is only singular at the boundaries of the integration region. Applying the following sequence 
of mappings 



A 2 i-> a(A 2 , A 3 , 1) 
A 4 a(A 4 , A 2 A 3 , 1) 
A 2 i-> a(A 2 , Ai, 1) 

leads to a factorized form. Numerically we then obtain 



h [J = 1] = A^4 



-201.16(3) 



68.426(4) 12.027(8) 1.0397(1) 0.75 



+ 



+ 



A last class of singular integrals with four denominators is: 
J(Pi,P2,P3,Pa) 



(6.9) 



(6.10) 



h 



\A-J 



S13S23S14S24 

J(Pl,P3,P2,P4) 



I 



S12S34S14S23 
d$ 4 (J(pi,P3,P2,P4) + J(P2,P4,P1,P3))- 



S24 



S34S12S23(S13«23 + S14S24) 

The term in parenthesis in the last denominator is given by 



(6.11) 



S13S23+S14S24 = A1A1A2 



2A2A3A4A4 + A 3 (Al + \\) + 2cos(7rA 5 )(2A 4 - l)V A2A3A3A4A4 

(6.12) 

We let d/3 = \2dI3 + \2dI3 and apply 



to \2dI3 and 



A 2 ^ a(A 2 , A 3 , 1) 
A 3i4 a(A 3j 4,A 2 ,l) 

A 4 i-> a(A 4 , A 3 , 1) 



(6.13) 



(6.14) 



to \2dlz- After these transformations, both integrals have only factorized singularities and 
can be evaluated numerically. We obtain: 



h [J = 1] = A^4 



-292.54(4) - 217 -° 3 °( 9 > - 52 - 76 5 2 ) + 6 - 931 ^ 3 ) + A 



(6.15) 



Other integrals with four denominators are mapped to ii,2,3[</] with simple relabeling 
of the partonic momenta. Integrals with fewer denominators are simpler and the non-linear 
mappings for them are obvious upon inspection, following the examples of Ref. [74]. One 
case which requires special attention arises when a denominator is raised to the second 
power in the squares of diagrams where a gluon splits into a qq or gg pair as has already 
been discussed in the literature (see, for example, [68,83]). These quadratic singularities 
are fake. A combination of terms, such as 



2m 2 H f S14 s 2 4 



'34 



S134 s 234 



(6.16) 
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is the one which emerges in the physical matrix elements. After a non-linear mapping 

A 2 ->a(A 2 ,A 3 ,l) (6.17) 

to factorize an overlapping singularity, the integrand is also free of quadratic singularities 
in the integration variables Aj. 



7. Numerical results 

In this section we present our numerical results for the decay rate ^n^blx VJ] f° r selected 
infrared safe observables J. We select the value of the renormalization scale fx = ran- 
First, we compute the inclusive decay width. Our numerical result is, 



-pNNLO _ r LO 
1 H->bb ~ 1 H^bb 



( OL s 
1+ — 
V TT 



5.6666(4) 



fCts_ 
V TT 



29.14(2) +0(a 3 s ) 



with 



,LO _ vtmHN 



L H^bb 4yr • 

This is in agreement with the known analytic expression [1] 



^NNLO 
H->bb 



r 



LO 
H^bb 



1 + 



TT J 



5.6666666.. + 



tty 



TT 



29.146714.. +0(a 3 s 



(7.1) 



(7.2) 



(7.3) 



With our numerical program, we can compute the decay rate for arbitrary infrared-safe 
observables. As an example we present our results for the 2,3 and 4 jet rates using the 
JADE algorithm [84] with parameter y cut = 0.01: 

(^) 2 94.1(l)+0(^) 



r^(4JetRate) = T_ 
r™- 6 (3JetRate) = Tf^ 



r™f (2JetRate) = T LO 



H^bb 



+ 
+ 
1 - 



(^) 19.258(4) + (^) 2 241(2) + 0{a\) 
13.591(6) - 307(2) +0{a 



TT 



(7.4) 



We have checked that our results for the jet-rates add up to the inclusive rate for general 
y cut values. 

In Fig. 1, we present the leading jet energy, E max , in the rest frame of the Higgs boson, 
for events with two jets (y cu t = 0.1). This is a new result which can only be obtained with 
a fully differential NNLO calculation and cannot be inferred from the inclusive decay width 
and NLO calculations. At leading order, E max = At higher orders, there is a range 
of jet energies which are kinematically allowed. We choose a value of a s {mz) = 0.118 at 
the Z boson mass and evolve consistently through LO, NLO and NNLO up to the Higgs 
boson mass, which we assume to be mn = 120GeV. 

The above numerical results demonstrate the applicability of our method to physical 
processes. A number of phenomenological studies which are relevant to the searches of the 
Higgs boson can be made. It is easy to interface our numerical code with a Monte-Carlo 
code for the production of a Higgs boson. We we will present complete phenomenological 
studies pursuing this direction in the future. 
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Figure 1: The energy spectrum of the leading jet in the decay H — » bb in the rest frame of the Higgs 
boson through NNLO. The jet clustering is performed with the JADE algorithm with y cut = 0.1. 

8. Conclusions 

In this paper, we present a first physical application of a new method for the factorization 
of overlapping singularities in phase space and loop integrations. We compute the fully 
differential decay width of a Higgs boson to a bottom-quark pair. We produce the required 
tree, one-loop and two-loop amplitudes with standard Feynman diagrammatic methods. 

Our article focuses on the phase space integrations which emerge at NNLO. We apply 
non-linear mappings to factorize all overlapping singularities in all real- virtual and double- 
real integrations. Consequently, we perform the expansion of all integrals in the dimension 
regulator e with simple subtractions. The formalism allows for the computation of the 
decay rate for arbitrary physical observables. 

We verify that we can reproduce the known results for the NNLO inclusive decay 
width and compute the differential two, three and four jet rates with the JADE algorithm. 
We also present the leading jet energy distributions, a new result that cannot be inferred 
from previous calculations. In the future, we will interface our NNLO Monte-Carlo to 
Monte-Carlo programs for the production of a Higgs boson. 

We believe that our method for NNLO computations is powerful and suitable for 
complicated physical processes. We are looking forward to further apply our method on 
other interesting processes at hadron colliders. 
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A. Double-real matrix elements 

The function B{pi,p2,Pz-,Pi) occurring in the H — > bbbb squared amplitude can be ex- 
pressed as 

t-j / , _i n _i . _i , - _i „ Su (S13 2 + 2S12 2 + 2si 2 Si 3 ) 
B(jPl,P2,P3,PA) = 2s 3 4 -2S134 - 4 5 234 + 2 s 124 -2 i ^ 

S14S34S124S234 

+ 4 fi 12 g 23 _ 2 512 (^23 + 2^13 + gig) 4 ^12^23 

S34S134S234 S14S34S124 S14S124-5134 

/ 4 5i 3 5i 4 + 2 g i3 2 + 4 S12S13 + 2 Su 2 + 2 Sl 2 Sl4 -2s 2 4 - 4si 3 

I S34S124S234 S14S34 

, -2 s 2 4 + 2 S12 - 2 S23 4g 14 - 4s 2 3 -4 s i2 - 4 s 14 - 4 s 13 



•Sl4«234 S124S134 S134-S234 

| 2fi 14 +4si 3 + 2fi23 2 glg£l3 -2s 2 4 + 2^12 _ 2 «13 (-«23 + S13) 



S34S124 Si4S 3 4Si24 S14S124 Si4S 3 4S234 

^ -4 513 - 2^14 + 2s 2 3 4 5l25l 3 + 2 524^23 + 2 524 2 -4£ig - 4 524 ~ 4 523 

•S 34 S234 S14S124S234 S134 2 

2s 2 4 + 2s 23 . -1,„ -512.S13 2 - -1,0 -1 , 2si2 + 2s 24 

H h 2 S14 1 + 2 8 S134 + 2 S234 H 

S34S134 S14S34S124S234 S14S134 

-2 si24 -1 ^ e 2 + ^ - 4 su" 1 + 6 S34" 1 + 6 si 34 _1 + 4 s 2 34 _1 - 2 s 12 a~ 1 

+2 gl2gl3 2 _ 4 512^23 _ 4 gl2^23 



S14S34S124S234 S34S134S234 S14S124S134 

_ 2 g 12 (-S23 ~ 2gj 3 + gig) 4si2 + 4g 2 4 + 4 S 2 3 
S14S34S124 -S134 2 
-6g 2 4 ~ 8gl 3 -6 gig ~ 2 5 23 + 6 5 34 + 8 S13 6 513 + 6 5 24 ~ 2 5 i4 + 2 S 2 3 



Si 4 S 3 4 S14S134 S34S134 

6513 +4 S 2 4 + 10 512 -2 5 2 3 ~ 4 5 24 + 2 5 i4 + 4 Sjg -4 5i 3 - 2 S U - As 12 



S14S234 S34S234 S134S234 

-2 S 2 4 - 2 S12 6 523 ~ 8 S12 - 8 514 + 6 $13 ~ 2 S13 - 4 S23 + 2 5 M 



^145124 S345124 S124S134 

4 514 - 4 513 - 8S12 , Sl3 («12 + S24 + S23) , ,513 (512 + 524 + 523) 



+ 4 " — ^+4 



S124S234 Sl4«134 2 S34S134 2 



4 S13S23 + 2 S12S23 - 2 S13 2 - 2 5 2 3 2 



S14S34S234 

2 Si2 2 + 4 S 2 3 2 + 8 S24S23 + 4 5i 2 523 + 6 5 2 4 2 + 8 S12S24 

S14S234S134 
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| -4 S12S23 + 6 5i4 + 4 S12S14 + 2 S23 + 2 S12 - 8 S23S14 2 5 2 4 + 2 S24S23 - 4 si 2 z 



S34S124S134 S14S124S234 

-4 S12 2 - 6 S12S14 + 2 S13 2 - 2 s M 2 \ 4 S12 + 4 s 23 + 8 ai 3 
H e H 

«34'Sl24S234 / s 14 s 34 

| -6 513 ~ 2si4 - 4s 2 4 2 g 2 3 ~ 2 5 2 4 + 4 gig - 2 S34 - 4 gig 

■534^134 514S134 
| -4 5 23 - 10 gig + 4 S 2 4 ~ 6 513 4fii2 + 4si 3 + 2s 2 3 + 2 514 

S14S234 S 34S234 
| 2s i4 +4si 3 +4s i2 2s i2 + 2s 2 4 -6 513 ~ 6 S23 - 4 512 2 513 ~ 2 5i4 + 4 5 2 3 

5l34«234 5145124 534S124 51245134 

| 4 5i4+4 5i 3 + 8 5i 2 d 5i 3 ($12 + ^24 + S23) { S 13 (gig + 5 24 + 5 2 3) 

51245234 514S134 2 S34S134 2 

| -2 5 2 3 2 ~ 2 5i 2 5 23 + 2 5i 3 2 - 4 5i 3 5 2 3 + 4 S i2 2 + 4 S 2 4 2 + 4 5i 2 5i 3 - 4 5 2 45l3 ~ 4 S12S24 



5145345234 

-2 Si2 2 - 2 S 2 4 2 - 4 5 2 3 2 - 4 S12S24 - 4 524523 - 4 S12S23 



51452345134 

-2s 23 2 - 2S12 2 + 4 5 23 514 - 2 5 M 2 , 4 Si 2 2 - 2 S24S23 ~ 2 5 24 2 



+ 



S3451245134 514S1245234 
| -4 5i 3 5l4 - 8 5l2 2 ~ 2 5l 3 2 - 2 5 M 2 - 8 5i 2 5i 3 - 6 5i 2 5i 4 ^ 

S34S124S234 



A 2L, S = 20 s 234 _1 + 20 s 134 ~ 1 + 32 s^ 1 + 24 524" 1 + 24 su" 1 + 32 S23" 1 

3 513 - 6 m H 2 + 3 S14 + 6 534 3 s 2 4 - 6 m# 2 + 3 s i4 + 3 s 34 

+4 h 5 

524523 S13S23 

+4 6 m H 2 ~ 3 513 + 3 534 + 3 «24 4 -524 ~ 514 + 5 34 + 4m g 2 

5234 514 5 2 34 S13 

4 —3 S14 — 6 rriH 2 + 3 S23 — 3 534 ^ 3 523 + 6 s 34 — 6 m# 2 + 3 s 2 4 

5134 524 5i 3 Si4 

-5s 34 - 8m H 2 - 4s 2 4 - 4 5i4 , m H 2 , _ mjf 2 

+0 h o j + o j 

5134 S234 5134 S234 

_ 4 8 m H 2 - 5 5i4 - 7 5 34 - 5 523 4 -5 2 4 ~ gl4 + 5 34 + 4m g 2 

524513 5134 5 2 3 

_ 4 - 5 5 2 4 + 8 m H 2 - 5 513 - 7 S34 _ 4 -5i4-2 534 + 5i3-3 m H 2 

514523 s 234 524 

+8 3 523 + 3 513 + 3 5 34 ~ 6 fflj? 2 _ ^ -3 TTT-ff 2 - 2 5 34 - 5 2 4 + 5 2 3 

524514 5134 S14 

^ ^ (m g 2 + 534) (-2m H 4 ~ 5 34 2 - 2m g 2 5 34 ) 

5234 5134 5i 3 S 2 3 

^ ^ (m g 2 + 534) (-2m g 4 - 5 34 2 - 2mjj 2 5 34 ) 

5234 5134 514524 

| 4 (-^34 + mg 2 ) (-2 m g 4 - 5 34 2 + 2 m H 2 5 34 ) 
S13514S23S24 
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_^3m H 2 s l3 - S34 2 - s l3 2 -4m H 4 - 3m fl 2 s 34 + si 3 s 34 



-4 
-4 
-4 
+4 
+4 
+4 
+4 
+4 
+4 
+4 



2 
-S34 


— 4 m# 4 - 


•S234S24S14 

- S24 2 + S24S34 + 3 iriH 2 S24 — 3 mH 2 ss4 




S14 — Si4 2 


S134 S13S23 
+ S14S34 — 3mn 2 S34 — 4m# 4 — S34 2 




— Attih - 


s 234 • s 13 s 23 

-3m H 2 S34 - S34 2 - s 2 3 2 + 3m H 2 s 2 3 


-S34 2 


-Am H 4 - 


S134 S14S24 

- S24 2 + S24S34 + 3 mH 2 S24 - 3 m H 2 S34 


-«24 2 


~ 3s 2 4«34 


«134 513^234 

— Auih 4 — 3s 3 4 2 — Qm,H 2 S34 — 3mH 2 - 



S234 «134 «14 

— 3S34 2 — 3^14534 — 4m/f 4 — 2>m,H 2 si4 — S14 2 — 6mjy 2 S34 

S134 «234 S24 

3m H 2 s u - S14 2 + 514^34 - 3m g 2 g 34 - 4m g 4 - s 34 2 
S234 S134 S23 

4 m H 4 + s 2 3 2 + 3 S34 2 ~ 3 m H 2 S23 - 6 m g 2 g 34 + 3 S23S34 

S24S13S14 

4 m H 4 + 3 S34 2 + s 2 4 2 - 6 m H 2 S34 + 3 S24S34 - 3 m H 2 S24 

S13S14S23 

-3 m H 2 si4 + 4 ran 4 + S14 2 + 3 S14S34 + 3 s 34 2 - 6 mH 2 «34 

S24S13S23 

-6 m H 2 S34 + 4 m H 4 + 3 S13S34 + S13 2 - 3 m H 2 si3 + 3 s 34 2 



•S24S14S23 



A Hbb 99 = -16s234- 1 - 32 Sl 3- 1 -24 Sl 4- 1 -16 Sl 34- 1 - 24 S 24- 1 - 32 S23 - 1 

, mfl 2 -, c m H 2 .Sl4 + 2s 34 ,S 2 4 + 2S 3 4 
-16 ^ - 16 5- - 4 4 

S234 S134 S234S24 S134S14 

+8 -3 S24 - 3 g 34 + 2 m H 2 - 3 s 14 4 -5 g 34 - 2 s 13 + m H 2 - 2 s 14 

Si 3 S 2 3 S24S23 

+4 -3 m H 2 - S34 + gi4 + g 2 4 4 - 3 g 2 4 - 3 5 34 - 3 m g 2 + 4 gig 
•S134S23 S234S14 
-3 m^ 2 + 8 S34 + 6 s 2 3 + 6 S14 3 m# 2 + 3 s 34 + 3 S14 - 4 s 23 

S24S13 S134S24 

+4 -3 m g 2 - s 34 + S14 + g 2 4 _ 4 -3 m g 2 + 8 S34 + 6 s 13 + 6 g 2 4 

S234S13 S14S23 

l0 5s 34 + Am H 2 + 4si4 + 4s 2 4 , -m H 2 + s u . -m H 2 + s 24 

+0 4 — 4 

S134S234 S234S23 s 134Sl3 

2 ran 2 - 3 S23 - 3 si 3 - 3 s 34 -2 s 24 - 5 s 34 - 2 s 23 + 2 

+0 h 4 

S24S14 S 13 S U 

s 13 2 - m H 2 s 13 - S13S34 + m H 2 S34 + s 34 2 
S234S24S14 
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S24 2 ~ m H 2 S2A ~ 524^34 + «34 2 + m H 2 S 34 
S134S13S23 

-^14^34 + m H 2 s 3 4 + s 34 2 - m g 2 Si4 + gj 4 2 
S234S13S23 

m H 2 S 3 4 + S23 2 ~ ^23^34 ~ m H 2 S2 3 + S 34 2 
Sl34«14 s 24 

| 4 -S14.S34 + rn H 2 S 3 4 + S34 2 ~ W-H 2 gl4 + S U 2 
S234S134S23 

| ^ s 24 2 ~ m H 2 S24 ~ S24S34 + S34 2 + m H 2 s 3 4 

S134S13S234 

| ^ s 24 2 + 3 s 34 2 + m H 2 s 2 4 + 3 S24S34 + 2 m H 2 S 3 4 
S234S134S14 

^ 4 2 m H 2 s 3 4 + m H 2 si4 + si4 2 + 3 s 34 2 + 3 S14S34 
S134S234S24 

| 4 ~gi3 2 + 2 m H 2 s 3 4 + m H 2 si 3 - 3 si 3 s 34 - 3 s 34 2 
S24S14S23 

| ^ -3 S23S34 + m H 2 s 23 - s 23 2 - 3 g 34 2 + 2 m H 2 su 
S24S13S14 

| ^ m H 2 S 2 4 + 2 m H 2 S 3 4 ~ 3 S24S34 - g24 2 - 3 s 34 2 
■Sl3«14S 2 3 

-3 s 34 2 - s i4 2 + m H 2 si4 + 2 m H 2 s 3 4 - 3 S14S34 



+4- 



S24S13S23 

2 



| 1 (-£34 + m g 2 ) g 3 4 2 | 1 (m g 2 + S34) S34 2 | i (m H 2 + S 34 ) g 34 2 ^ _ 
S13S14523S24 S234S134S13523 S234S134S14S24 



A mb 99 = " 4S 23" 1 - 4S13" 1 - 8S134" 1 - 8S234" 1 



l0 m H 2 o m H 2 2m H 2 + 2s 3 4 s 34 + m H 2 - s 23 

+0 x + o 5- + 8 4 

S134 «234 -Si 34 S2 3 4 S134SU 

_ 4 ^34 - S13 + m H 2 4 ~2g 34 + S13 - S 2 4 _ 4 3 g 34 + S i4 + g 2 3 

S234S24 S234S14 S24Sl 3 

+4 £14 + S24 ~ S34 _ 4 £14 + 2 5 34 ~ S23 4 $14 + ^24 ~ «34 

Sl34«23 S134S24 S234-S13 

_ 4 S24 + 3 S34 + £13 _ m g 2 - £14 _ m g 2 - S 2 4 

514^23 S234S23 S134S13 

+4 -S13 + ^g 2 ~ 3 s 3 4 - su 4 -g 2 4 ~ 3 s 34 + m H 2 - s 23 

S2AS 23 S 13 Si4 
_ 4 m H 2 S 3 4 - S 13 S 3 4 + ^34 2 _ 4 «34 2 ~ ^24^34 + m H 2 S 34 

S234S24S14 Si 34 Si 3 S 2 3 
_ 4 S34 2 + m H 2 S 3 4 - SuS 3 4 _ 4 -^23^34 + m H 2 S 3 4 + S34 2 

S234S13S23 Si 34 Si 4 S24 
+4 S34 2 ~ ^24^34 + m H 2 S 3 4 4 524^34 + m H 2 S 3 4 + 2 S 34 2 

S134S13S234 S234S134S14 
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+4 S14S34 + m H 2 su + 2 s 34 2 ^ g 34 2 + m H 2 S34 - s u s 34l 

S134S234S24 -S234S134S23 

+4 -2 s 34 2 - §23^34 + rn H 2 s 34: 4 -^24^34 + m H 2 s M - 2 g 34 2 

S24S13S14 S13S14S23 
+4 -2 s 34 2 - si 4 s 34 + m H 2 s 3i 4 -2 S34 2 - gi 3 s 34 + m H 2 s u 

•S24S13S23 S24S14S23 

, (-S34 + m-H 2 ) S34 2 , (w-h 2 + S34) S34 2 , [mu 2 + S34) s 34 2 ,. 
+4- |-4- |-4- (A.4) 

Sl3 s 14S23S24 S234Sl34«13'S23 52345134514^24 



Klhgg = 8 Sl4_1 + 16 S2 ^ + 8 S24 ~ 1 + 16 + 40 Sl34_1 + 40 S234 " 1 

m H 2 (-2 m H 2 + s 34 + su) m fl 2 (-2 m fl 2 + s 34 + S24) 

+0 h 8 

S234524 s 13 S134S14S23 

m^ 2 (2 mf/ 2 - si 3 + S24) _ "iff 2 (2 m// 2 + s u - s 23 ) 
— 8 — 8 

S234S14S23 S134S13S24 
niH 6 itih 6 „si4 + 2m// 2 m# 2 



+16 + 16 + 8 — — + : 



S134S13S234S24 S234S134S14523 S234S13 5134S13 

-3 m H 2 + S23 - S24 - S34 2 S13 + s 34 + 2 S24 - m H 2 

+ ' 



S134S14 S14S23 

,i3-. 3 4-3m H 2 -,i4 , - m H 2 m H 2 s 2A m H 2 s u 

—o h o o k + o 5 

S234S24 5234^23 S23S234 «134 «14 

^"^34 m H 2 s u , n -V-4m^ + 3mA4 

+ 8 2~ 8 ~ 2 + 



+8 



«234 S24 S134 S13 -S134S13S234 

-2^24^34 ~ 3m g 2 s 34 - 4m g 4 - s 2 4 2 - S34 2 ~ 3m H 2 s 2 4 
•S234S134S14 

-3 m H 2 s 3 A - si4 2 - 2 S14S34 - S34 2 - 4 m H 4 - 3 m H 2 su 



•S134S234S24 



Sm H Z Su - 4 TOtf 4 - Si4 Z 0/1 "T,// 2 , 0/l «7,tf 2 
+0 h Z4 » + Z4 k 

S234S134S23 -5134 S 2 34 

-S34 + si 3 - 3 m H 2 - 2 s 2 4 -"iff 2 + S34 + 2 S14 + 2 s 23 

+ 1 



+8 



S234S14 S24S13 
s 2 4 + 2m H 2 -S34 - 2 gi4 - 3 m H 2 + s 2 3 
S134S23 S134S24 



, . r -6 m// 2 - 4 S24 - 4 S14 - 3 s 34 

+16 (A. 5) 

S134S234 



= "8 su" 1 - 16 s 2 3 _1 - 8 S24 -1 - 16 sis" 1 - 48 s^ 1 - 48 S234" 1 

io m H 2 (-2s 34 - su) m H 2 ( — 2 S34 — S24) m H 2 (s 13 - 2 s 24 ) 

+0 h o o 

S234S24S13 S134S14S23 S234S14S23 

mij 2 (-2 S14 + S23) S13 o -m H 2 + S34 - «23 
— o h o o 

■S134S13S24 S234S23 -S134S14 

n -S13 - m// 2 + S34 , , r . m H 2 s 2 4 , . m H 2 s 3 4 



+ 16 , - 16 



S234S24 S23S234 2 Si 34 2 Si4 
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, R m H 2 s M m H 2 su „„ V QO m H 2 

— lb k h lb 7, 62 » — 62 » 

«234 «24 Si34^Si 3 S 1M Z S 23 4 

+8 S24 2 - m H 2 s 2 4 8 2 S24S34 + g34 2 + S24 2 + m H 2 s 2 4 + m H 2 s 3A 

•Sl34Sl3 s 234 S234«134Sl4 

S34 2 + 2 S14S34 + m H 2 s 1A + S14 2 + m H 2 s M s 14 2 - m H 2 s 14: 

+5 h o 

s 134 s 234-524 S234S134S23 

-S13 + m^ 2 + 2 s 2 4 + S34 , _ -2 S14 + m^ 2 - s 34 - 2 s 23 

— o h o 

•S234S14 S24S13 

, Q -2s24- 3 m H 2 , _ «23 mtf 2 - S 23 + S34 + 2 S14 

+0 h o o 

S134S23 S134-S13 Si34S 2 4 

+8 -3m H 2 - 2su 8 -g 3 4 - 2 g 24 - 2 s 13 + m H 2 

S234S13 S14S23 
, , r . 4s M + 4s 2 4 + 2m H 2 + 3s 34 

+lb (A.b) 

S134S234 



R 2) S34 + S 24 _i _i _i _i 

B Mhhnn = ~ 8 8fi 23 -8si 3 -I6S134 -lbs 23 4 

m H 2 s 34: . rn H 2 s 3 4 m H 2 S24 m H 2 s 14: 
+0 r o o o ■ 



S234S24S13 -S134S14S23 S234S14S23 -S134S13S24 

-S13 + m// 2 S34 + S14 m H 2 S24 

+0 o o ;r 

S234S23 S234S24 S23S234 

m H 2 su . Q m H 2 s u m H 2 s 14 

+0 „ h o k o „ 

S134 «14 «234 S24 «134 «13 

. o -S24S34 + m H 2 s u s 34 2 + s 24 s 3 4 + m H 2 s u 

+0 h o 

•Sl34Sl3'S234 ^234^134514 

, c s 34 2 + m H 2 s 34: + S14S34 -S14S34 + m H 2 s 34: 

+0 h o 

S134S234S24 S234S134S23 

m// 2 , m// 2 s 34 + m# 2 + s 2 4 - s 13 

+0 ^ + o k — o 

S134 s 2 34 S234S14 

. -S23 - S14 - S34 + "I// 2 . o S 24 + S14 , „ -S23 + m H 2 
+0 h o h o 

S24S13 Sl34«23 Sl34«13 

m# 2 + S34 - S23 + S14 , s 24 + Sl4 
— o h o 

S134S24 S234S13 

m fl 2 - s 34 - S13 - s 24 , , „ 2 mf/ 2 + 2 s 34 , . 

+8 h 16 (A. 7) 

S14S23 S134S234 



^(0) _ g 2% + 4m g 2 + gi4 - S13 _ 4 -£14 - 2 S13 + 4 m g 2 - 4 s 2 3 



S234S34 S24S34 

_ 4 -2 S14 - 2 S23 + 2 m g 2 - s 2 4 _ 4 2 m g 2 - 2 gig - g i4 - 2 s 2 4 

■S13S34 S23S34 

, , a -m H 2 s 2 A ~ 4 m^ 4 - 2 s 2 4«i4 - ^u 2 - «24 2 - m H 2 s M 

+16 

S134S34S234 



- 20 - 



-S24 2 + 2 m_ff 2 s 2 4 - 2 m// 4 - si 4 2 + 2 mj^su - 2 S24S14 

S234S34S13 

4 2 g 2 4.si3 ~ 2 m g 4 - s 2 4 2 - si 3 2 + 2 m H 2 si 3 - 2 m H 2 s 2 4 

S234S34S14 

_ 32 g24Sl4^g 2 _ 4 ~g 2 4 ~ 2 S 2 3 + 4 m g 2 - 4 S13 



-4 



Sl34«34 z S234 S14S34 

2 ^14^23 ~ g23 2 + 2 m H 2 s 2 3 - 2 m g 4 - g i4 2 - 2 m H 2 s 14: 
S134S34S24 

-S24 - 4mn 2 - 2S14 + S23 , 1n m H 2 S24: 



+ 16 ■ 

■S134S34 S234 S34 

2 m H 2 S2i + 2 m H 2 si3 - 2 m g 4 - s 2 4 2 - 2 S24S13 ~ si3 2 
S14S23S34 

, -, a S2A 2 m H 2 1c ,m H 2 s u , i„s 14 2 m H 2 
+16 — ^ ^ + 16 ^ h 16 — ^ tt 

S34 S234 «134 S34 S34 S134 

-2 514^23 + 2 m H 2 S23 - 2 mjf 4 - s 2 3 2 + 2 m H 2 su - s\ 2 

S24S13S34 

_ 4 -S24 2 + 2 m H 2 S2A - 2 m g 4 - s u 2 + 2 m H 2 s 14: - 2 g 2 4^i4 ^ ^ 

S134S34S23 



£,(1) _ g -S14 - 2fi24 + g!3 ^ «24 + 2 g 2 3 + 514 ^ Sl4 + S 2 4 + 2 S 13 
Hbl " S234S34 S13S34 ' S23S34 

S 2 4S14"Z J H' 2 , 2s 23 +4S13 + S 24 , ~3 514^23 + «23 2 + «14 2 



+32 4 — ^ ^ - 4 

S134S34 S234 S14S34 S134S34S24 

2 S14 + s 24 - S 2 3 -. a rn H 2 s 2 4 . 3 S 24 Sl3 + S13 2 + s 24 2 
—8 16 5 4 

S134S34 «234 S34 -S14S23S34 

, r s 2 A 2 m H 2 . 4«23 + «14 + 2si 3 -.„rn H 2 su 
-16 — k k - 4 16 5 

S34 S234 S24S34 S134 S34 

m H 2 S 14: + 2 S 2 4«14 + rn H 2 S2A + S14 2 + «24 2 , e si4 2 "iH 2 

+ 16 16 K 7J 

S134S34S234 S 3 4 Z Si34 Z 
_ 4 S 2 4 2 + «14 2 + 3 S24«14 ^ S24 2 + Sl4 2 + 3 S24«14 
•S234S34S13 S134S34S23 
4 g i4 2 + 3 S14S23 + «23 2 4 Sl3 2 + S 2 4 2 ~ 3 S24-S13 ^ Ag ^ 
524S13S34 S234S34S14 
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